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Observation of muon intensity variations by season with the MINOS far detector 
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The temperature of the upper atmosphere affects the height of primary cosmic ray interactions and the pro- 
duction of high-energy cosmic ray muons which can be detected deep underground. The MINOS far detector 
at Soudan MN, USA, has collected over 67 million cosmic ray induced muons. The underground muon rate 
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measured over a period of five years exhibits a 4% peak-to-peak seasonal variation wliich is higiily correlated 
with the temperature in the upper atmosphere. The coefficient, qt, relating changes in the muon rate to changes 
in atmospheric temperature was found to be: ar = 0.873 ± 0.009 (stat.) ±0.010 (syst.). Pions and kaons 
in the primary hadronic interactions of cosmic rays in the atmosphere contribute differently to ar due to the 
different masses and lifetimes. This allows the measured value of ay to be interpreted as a measurement of the 
K/tt ratio for Ep >7 TeV of 0.12lo Qg, consistent with the expectation from collider experiments. 

PACS numbers: 95.55.Vj, 98.70. Vc, 98.70.Sa, 14.44.Aq 



I. INTRODUCTION 

When very high energy cosmic rays interact in the strato- 
sphere, mesons are produced in the primary hadronic interac- 
tion. These mesons either interact and produce lower energy 
hadronic cascades, or decay into high energy muons which 
can be observed deep underground. While the temperature of 
the troposphere varies considerably within the day, the tem- 
perature of the stratosphere remains nearly constant, usually 
changing on the timescale of seasons (with the exception of 
the occasional Sudden Stratospheric Warming [1]). An in- 
crease in temperature of the stratosphere causes a decrease in 
density. This reduces the chance of meson interaction, result- 
ing in a larger fraction decaying to produce muons. This re- 
sults in a higher muon rate observed deep underground |l2|-|4j]. 
The majority of muons detected in the MINOS far detector 
are produced in the decay of pions, although the decays of 
kaons must be considered for a more complete description of 
the flux |5]. 

MINOS is a long baseline neutrino oscillation experi- 
ment |6, 7], with a neutrino source and near detector at Fermi 
National Accelerator Laboratory in Batavia, IL, USA, and a 
far detector at the Soudan Underground Mine State Park in 
northern Minnesota, USA. This paper describes cosmic ray 
data taken in the far detector, a scintillator and steel track- 
ing calorimeter located 0.72 km underground (2080 mwe, me- 
ters water equivalent) js*]. It has a 5.4kton mass and a 
6.91 X 10^ cm^sr |@] acceptance. Because of its depth, MI- 
NOS detects cosmic -ray muons with energy at the surface, 
>0.73 TeV. These high energy muons are mostly the result 
of the decays of the mesons produced in the primary hadronic 
interaction. This, coupled with the large acceptance, makes it 
possible to detect small seasonal temperature fluctuations in 
the upper atmosphere. The far detector is the deepest under- 
ground detector with a magnetic field, allowing the separation 
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of particles by charge. 

The MINOS data are correlated with atmospheric tem- 
perature measurements at the Soudan site provided by 
the European Centre for Medium-Range Weather Forecasts 
(ECMWF) [10]. This temperature data set has higher pre- 
cision than any other used for the seasonal variation analy- 
sis iB l4l [lll - llTIl . The 67.32 million muon events used in this 
analysis were collected over five years, from August 1, 2003 
to July 31, 2008, a period that includes five complete an- 
nual cycles. The seasonal variations in muon intensity were 
compared to a theoretical model which extends the pion-only 
model of [3] to include the contribution from kaons. 



II. THE EXPERIMENTAL EFFECTIVE TEMPERATURE 
COEFFICIENT 

A. Experimental Intensity 

The underground muon intensity depends on the threshold 
energy Eth and the cosine of the zenith angle 9. The change in 
underground muon intensity variations as a function of tem- 
perature was derived following the formalism of |2, 3]. The 
change in the surface muon intensity, A/^(£', cos 6*) occur- 
ring at the MINOS far detector site can be written as: 

/•OO 

A/^= / dXW{X)AT{X) (1) 
Jo 

where AT{X) is the change in atmospheric temperature at 
atmospheric depth X, and the weight W{X) reflects the tem- 
perature dependence of the production of mesons in the atmo- 
sphere and their decay into muons that can be observed in the 
far detector A temperature coefficient a{X) can be defined 
as: 

a{X) = ^W{X), (2) 

where is the muon intensity evaluated at a given value 
of atmospheric temperature Tq. The phenomenological re- 
lationship between the atmospheric temperature fluctuations 
and muon intensity variations can now be written as: 

The atmosphere consists of many levels that vary continu- 
ously in both temperature and pressure. To simplify calcula- 
tions, the atmosphere is approximated by an isothermal body 
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with an effective temperature, Tcff, obtained from a weighted 
average over the atmospheric depth: 



Tcff — 



J^^dXT{X)W{X) 



(4) 



An "effective temperature coefficient", ar can then be defined 



off 



I?, ./o 



dXW{X). 



(5) 



With these definitions in place, the relationship between atmo- 
spheric temperature fluctuations and muon intensity variations 
can now be written as: 



P 



ax- 



Trff 



(6) 



The configuration and geometric acceptance of the far detec- 
tor remain constant over time. Therefore, the rate, i?^ of 
muons observed in the detector is proportional to the incident 
muon intensity and varies with the effective atmospheric tem- 
perature as follows: 



ax 



(Tcff) ■ 



(7) 



In practice, the observed muon rates and the temperature data 
are averaged over the period of a day. The effective tempera- 
ture is obtained from a weighted average of temperature mea- 
surements obtained at a set of discrete pressure levels. 

The weight W{X) can be written as the sum + W^, 
representing the contribution of pions and kaons to the over- 
all variation in muon intensity. The weights W^'^ are given 



byfllSllia]: 



W'^{X) 



where 



K 



7 + (7 + 1) Si j,K{X) {{Eti, cos 0) /e^.Kf ' 

(8) 



K{X) 



(l-e-^/^--)A;,^/X- 



(9) 



The parameters ^ include the amount of inclusive me- 
son production in the forward fragmentation region, masses 
of mesons and muons, and muon spectral index; the input 
values are A^. = 1 and A]^ = 0.38 ■ r^/^r l>18t ,19,1 . where 
is the K/tt ratio. The parameters ^ reflect the rel- 
ative atmospheric attenuation of mesons; The threshold en- 
ergy, Eth, is the energy required for a muon to survive to a 
particular depth; The attenuation lengths for the cosmic ray 
primaries, pions and kaons are Ajv, Ajr and Ak respectively 
with 1/A^ ^ = I/Aat — I/Att^k. The muon spectral index 
is given by 7. The meson critical energy, e^^.K, is the meson 
energy for which decay and interaction have an equal proba- 
bility. Since the distribution has a long tail (Fig.|3]l, the value 
of (i?th cos 6) used here is the median. The values for these 
parameters can be found in Table |T] 



TABLE I: Input W(X) parameter values. 



Parameter 



Value 



A], 

rK/n 

Bk 
An 

Avr 

Ak 

{Eth cos t 

7 

en 

ex 



1 U8,i9J 

0.38 • r^7, [18, 19] 
0.149 [20J ±0.06 [2r] 
1.460± 0.007 |i8,i£] 
1.740 ± 0.028 [IS, 19} 
120g/cm2 [20] 
180g/cm2 [20] 
160g/cm^ [20] 
0.785±O14TeV 
1.7±01 [5] 

0.114±O003TeV [18, 19] 
0.851±O014TeV [18, 19] 



Since the temperature is measured at discrete levels, the in- 
tegral is represented by a sum over the atmospheric levels X„: 



EloAx„r(x„)(VK- 
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(10) 



where W^'^ is W^'^ evaluated at X„. The temperature and 
pressure vary continuously through the atmosphere. Fig. [T] 
(solid line) shows the average temperature from 2003-2008 
above Soudan as a function of pressure level in the atmo- 
sphere [10]. The height axis on the right represents the aver- 
age log-pressure height, the height of a pressure level relative 
to the surface pressure, corresponding to the average tempera- 
tures plotted here. The dashed line is the weight as a function 
of pressure level W{X), obtained from Eq.[8]and normalized 
to one, used to calculate the effective temperature. The dashed 
weight curve in Fig.[T]shows that the temperature fluctuations 
higher in the atmosphere have a greater effect on the produc- 
tion of muons visible at a depth of 2100 mwe. High energy 
mesons produced at the top of the atmosphere are more likely 
to decay, producing muons visible to MINOS, than those pro- 
duced lower in the atmosphere. Note that the expression used 
to calculate Tcff in the pion scaling limit, ignoring the kaon 
contribution, is the same as the MACRO calculation f^. The 
effective temperature coefficient, ax, is a function of both the 
muon threshold energy and the K/tt ratio. As the energy in- 
creases, the muon intensity becomes more dependent on the 
meson critical energy, which in turn is proportional to the at- 
mospheric temperature. The effective temperature coefficient 
thus reflects the fraction of mesons that are sensitive to atmo- 
spheric temperature variations, and for energies much greater 
than the critical energy, the value of ax approaches unity. At 
the depth of the MINOS far detector, the vertical muon thresh- 
old energy lies between the pion and kaon critical energies. 
Therefore, because the muon energy is close to the parent me- 
son's energy, a larger K/tt ratio results in a smaller value of 
ax- 
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FIG. 1: The five year average temperature at various pressure 
levels (solid line). The range is from 1000 hPa 
(lhPa= 1.019 g/cm^), near Earth's surface, to 1 hPa (nearly 
50 km), near the top of the stratosphere. The height axis on the 
right represents the average log-pressure height corresponding 
to the average temperatures plotted here. The dashed line is 
the weight as a function of pressure level (X) used to find Tcb. 
The weights are determined by Eq.jS] normalized to one. 




Ays) 



FIG. 2: The time between consecutive cosmic ray muon 
arrivals, fit with a Poisson distribution. The fit gives 
X^/NdoF = 55.2/68; 0.4692 ±0.0001 Hz (from 

slope). The Poissonian nature of the muon arrival times 
demonstrates the absence of short timescale systematic effects 
on the data. 



B. The Data 

The muon data for this analysis were accumulated over a 
five year span, beginning on August 1, 2003. Data quality 
cuts were performed to ensure a clean sample of muons (Pre- 
Analysis cuts) JKl 

1 . Require that all detector readout and sub-systems were 
functioning normally 



2. Remove runs with anomalous cosmic ray rates, greater 
than 1 Hz 

3. Remove events that had many hits assigned to incorrect 
channels (properly demultiplexed [5]) 

4. Remove muons induced by NuMI beam interactions 
with timing cuts |@]. 

After all cuts were applied the initial sample of 68.66 mil- 
lion muons was reduced to 67.32 million muons |18]. A 
plot of the time between consecutive muon arrivals in the 
MINOS data is shown in Fig. |2] The distribution is well 
described by a Poisson distribution ifTsl l23ll with mean rate 
(Rf,) = 0.4692 ± 0.0001 Hz, demonstrating the absence of 
short-timescale systematic effects on the data. The average 
muon rate was calculated for each day by dividing the num- 
ber of observed muons by the detector livetime. 

The energy spectra for the observed muons can be seen in 
Fig. [3] The solid line is Eth cos 6, which was used to deter- 




2000 4000 6000 8000 10000 

E(GeV) 

FIG. 3: A plot of the energy spectra observed in the far 
detector The dashed line is Eth cos 0, which was used to 
determine the value used in Eq.[8] The solid line is the 
distribution of muon surface energies, in far detector. Also 
shown are Eth cos 6(CS) (dot-dash line), the distribution of 
Eth cos 6 after charge-separation cuts have been applied, and 
E^{CS) (dotted line), the distribution of E^ after the 
charge- separation cuts (see Sec. Ill Cb have been applied. Note 
that the charge-separation cuts have been applied, but the 
distributions shown include both muon species. 



mine the value used in Eq. [8] The dashed line is the distri- 
bution of muon surface energies, i?^, in far detector, which 
has a much longer tail than the distribution of threshold ener- 
gies. Also shown are Eth cos9{CS) (dot-dash line), the dis- 
tribution of Eth cos 9 after charge-separation cuts have been 
applied, and Ef^{CS) (dotted line), the distribution of i?^ af- 
ter the charge-separation cuts (see Sec. Ill CI) have been ap- 
plied. Note that the charge-separation cuts have been ap- 
phed, but the distributions shown include both muon species. 
The Eth cos 9 distribution is peaked, with a median value 
(£:thCOs6l)=0.795±0.14TeV. This distribution, with its rapid 
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FIG. 4: The daily deviation from tiie mean rate of cosmic ray 
muon arrivals from 8/03-8/08, shown here with statistical error 
bars. The periodic fluctuations have the expected maxima in 
August, minima in February. The hatched region indicates the 
period of time when the detector ran with the magnetic field 
reversed from the normal configuration. 
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FIG. 5: The daily deviation from the mean effective 
temperature over a period of five years, beginning when the far 
detector was complete, 08/03-08/08. The hatched region 
indicates the period of time when the detector ran with the 
magnetic field reversed from the normal configuration. 



TABLE II: Systematic errors on the experimental parameter 
inputs to ot . 



fall-off, reflects the rock overburden surrounding the far de- 
tector. 

The temperature data for the Soudan site was obtained from 
ECMWF, which collates a number of different types of ob- 
servations (e.g. surface, satellite and upper air sounding) at 
many locations around the globe, and uses a global atmo- 
spheric model to interpolate to a particular location. For this 
analysis, the ECMWF model produced atmospheric tempera- 
tures at 21 discrete pressure levels: 1000, 925, 850, 700, 500, 
400, 300, 250, 200, 150, 100, 70, 50, 30, 20, 10, 7, 5, 3, 2 
and 1 hPa (1 hPa = 1 .019 g/cm^), at four times, 0000 h, 0600 h, 
1200h and 1800h each day. The effective temperature. Toff, 
was calculated four times each day using Eq. [10] A mean 
value (Tcff ) and error was obtained from these four daily mea- 
surements. The ECMWF temperature data was cross-checked 
using the Integrated Global Radiosonde Archive (IGRA) of 
temperature measurements ll24ll . The distribution of the dif- 
ferences between ECMWF and IGRA temperature values at 
International Falls, MN was well described by a Gaussian dis- 
tribution with a = 0.31 K. 

Fig. m shows the percentage deviation in the mean daily 
muon rate, Ai?^, over the entire set of data, with statistical er- 
ror bars. A typical day at (i?^) = 0.4692 Hz yields -40,000 
muons, resulting in error bars of order 0.5%. The variation 
with season can be seen, with maxima in August and minima 
in February. These maxima peak at rates that are within 0.5% 
of each other. For the five year period (Tcff) = 221.93 K. 
The distribution of AToff over the data period can be seen 
in Fig. |5] with strong periodic seasonal correlation with the 
data. There is also striking correspondence between Fig. |4] 
and Fig.|5]for small term maxima and minima over a few days' 
span. 

A plot of Ai?^/ (Rf,) (ATcff) was produced (Fig. ^ for 
each day's Ai?^ and AToff data to quantify the daily correla- 
tion between rate and temperature. To find the value for ax, 
a linear regression was performed using the MINUIT Ii25il fit- 
ting package. This package performs a linear regression ac- 
counting for error bars on both the x and y axis using a nu- 
merical minimization method. The result of this fit is a slope 
of ax = 0.873 ± 0.009 (statistical errors only), and the cor- 
relation coefficient (R-value) between these two distributions 
is 0.90. 

The effects of systematic uncertainties were evaluated by 
modifying parameters and recalculating ax- Table HIl shows 
the difference in calculated ax for the modified parame- 
ters. The largest systematic errors are: a) the ± 0.06 un- 
certainty in meson production ratio f^P]; b) the ±0.31 K un- 
certainty in mean effective temperature, estimated by com- 
paring ECMWF temperatures at International Falls, MN, to 
those of the IGRA [24] measurements; c) the ±0.14TeV un- 
certainty in muon threshold energy, estimated from uncertain- 
ties in the rock overburden above the far detector. To esti- 
mate this uncertainty, the rock map was adjusted up by 10% 
and {Eth cos 9) was calculated, then the rock map was ad- 
justed down by 10% and (£'thCos6') was again recalculated. 
The uncertainty was then calculated from the difference be- 
tween {Eth cos 9) and these adjusted values. These system- 
atic errors were added in quadrature and are included with the 



Parameter Aqt 

meson production ratio, r^/Tr = 0. 149±0.06 [21J 0.007 

mean effective temperature, (Tcff )= 221.93±0.32 K 0.0051 

threshold energy, {Eth cos 6l)=0.795±0. 14 TeV 0.0048 

kaon constant, B],- = 1.740 ± 0.028 0.00046 

pion constant, Bl = 1.460 ± 0.007 0.000063 

Total 0.010 



error from the linear fit to obtain the experimental value of 
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III. DISCUSSION 
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FIG. 6: A plot of AR^/ (ii^,) as a function of AT^b / (Teff ) 
for single muons, fit by a line with the y-intercept fixed at 0. 
The fit has a x^/NdoF = 1959/1797, and the slope is 
QT = 0.873 ± 0.009. 



ar = 0.873 ± 0.009 (stat.) ± 0.010 (syst.). 



C. Charge Separated 



To obtain a sample of events with well-measured charge 
sign, further selection requirements were applied to the length 
and radius of curvature of muon tracks. These cuts, taken 
from previous investigations of the muon charge ratio at MI- 
NOS ISft, have the effect of reducing the energy distribution at 
Earth's surface of the selected muon sample. 

In all, 5.7% of the data set survived the cuts for both the 
forward and reverse field detector configurations. For the 
charge-separated samples linear regressions yielded effective 
temperature coefficients, arif^'^) = 0.79 ± 0.05 and 
axifJ--) = 0.77 ± 0.06 with x^/^Vdoj. of 1933/1758 and 
1688/1751 respectively. These numbers are consistent with 
each other, so there is no measurable difference between the 
temperature effect on /i+ and /i^. The value of the charge- 
separated ar is expected to be smaller than the previous ax 
with no charge separation because the selection cuts change 
the energy distribution over which the integration is performed 
to calculate ay. This can be seen in Fig. [3] with the most dra- 
matic difference between the all muon and charge separated 
distributions of E^cosO. This difference could produce the 
systematic offset observed between these values, and is dis- 
cussed further in the next section. 



A. Predicted ar 

The theoretical prediction of ar can be written as 101: 



1 



(11) 



Using the differential muon intensity ||20|] . 



dE„, 



/ V^iE, X)dX ~ A X 

"'0 

1 , 0.38r;^/, 



(I?) 



,1 + I.IS^ cos6'/e^ 1 + lAE^cosO/ex , 

and the MACRO approximation for the muon intensity iH, 
the prediction for ar can be calculated: 



aT 



1 l/eK + A],{DjDK?/e^ 



where 



Dtt,K — 



7 



7 + 1 l.li;thCos( 



1, 



(13) 



(14) 



Note that this can be reduced to MACRO'S previously pub- 
lished expression [ar)^ 131], by setting — (no kaon 
contribution). — 0.38 • rx/-K is the same as in SecUll 

A numerical integration using a Monte Carlo method was 
performed to find the predicted value of the seasonal effect 
coefficient, {ax)^, for the far detector A set of muons was 
generated by drawing values of E^ and cos separately from 
the differential intensity of muons at the surface, calculated 
in [20]. A random azimuthal angle, 0, was assigned to each 
event and combined with cos 6* and the Soudan rock over- 
burden map [5] to find the slant depth, S'(cos 6*, (/>), of the 
event. This was converted into the coiTesponding threshold 
energy, i^th, required for a muon on the surface to propagate 
to the far detector Events satisfying E^ > ii'th were retained, 
and the mean value of ax was found for a sample of 10,000 
events, giving (aT)p = 0.864 ± 0.024 for MINOS. When 
this calculation was performed using the lower energy charge- 
separated energy spectrum, the result is an (aT)p value that is 
lower by 0.015. This is most clearly seen in Eq.[T3j which is 
dominated by the leading 1/Z?^ term. As EthCosO increases, 
Z^TT goes to one. Any selection that reduces the Et\yCos9 dis- 
tribution will then reduce the expected ay- 

The systematic uncertainty on (aT)p was found by modi- 
fying the input parameters and recalculating ut- The dom- 
inant contributions were from: a) the ± 0.06 uncertainty in 
meson production ratio; b) the ± 10% in rock map uncer- 
tainty'; c) the ±0.1 uncertainty in muon spectral index; d) 



' The rock map is not a determination of the slant depth by geophysical 
means. It was created by measuring the muon flux coming from a partic- 
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TABLE III: Systematic errors on the theoretical parameter 



inputs to qt- 

Parameter Aqt 

meson production ratio, K/tt = 0. 149±0.06 [21] 0.020 

rock map uncertainty ±10% 0.013 

muon spectral index, 7= 1.7 ± 0.1 0.0031 

kaon critical energy, e/,- =0.85 1± 0.014 TeV 0.0014 

pion critical energy, e^=0. 1 14±0.003 TeV 0.0002 

Theoretical Total 0.024 



imental data points are taken from the MACRO survey and 
were calculated using a definition which excluded the contri- 
butions from kaons and were limited by temperature measure- 
ments up to 20g/cm^; when the MINOS result is recalculated 
with this definition the effective temperature coefficient de- 
creases to ax = 0.835. To compare the experimental values 
with the theoretical model, Eq.[T3] the expected effective tem- 
perature coefficient as a function of depth was calculated us- 
ing the numerical integration method outlined earlier, using 
standard rock and a flat overburden, and is shown in Fig. [T] 
as the solid line. There is qualitative agreement between the 
prediction and the experimentally measured values, but quan- 
titative comparisons would require recalculating the experi- 
mental values using the kaon-inclusive definition of effective 
temperature. The two dashed lines in Fig.|2]show the effective 
temperature dependence for the extreme pion-only and kaon- 
only predictions. Fig.|2]is illustrative only, as the dependence 
of the experimentally measured effective temperature coeffi- 
cient on the input K/tt ratio is not explicitly shown. 



B. Measurement of Atmospheric K/tt Ratio 

The uncertainty on the atmospheric K/tt ratio in the current 
cosmic ray flux models is of order 40% [2 1.1 . There has not 
been a measurement of this ratio with cosmic rays. Previous 
measurements have been made at accelerators for p-t-p colli- 
sions m, AuH-Au coflisions [29], Pb-nP coflisions SlU 
and pH-p collisions ll33ll . Many other older measurements are 
summarized in [32^. The experimental and theoretical values 
of ar can be combined to give a new measurement of the K/tt 
ratio for the reaction p + Aatm, with Ep >7 TeV. The thresh- 
old muon surface energy, £'th=0.73 TeV and the median muon 
surface energy, (£^^), is 1.15 TeV. On average, the muon en- 
ergy is one tenth the energy of its parent primary. The theo- 
retical ar depends directly on the K/vr ratio, as a consequence 
of the different interaction and decay properties of kaons and 
pions in the atmosphere. Since kaons and pions have different 
critical energies and attenuation lengths, the effective temper- 
ature also depends on the K/tt ratio, and therefore the exper- 
imental aT is a weak function of the K/tt ratio. By plotting 
the experimental and theoretical values of ot as functions of 
the K/tt ratio and finding the intersection of the two curves, a 
measurement of the K/tt ratio can be obtained. 

Fig.|8]shows the experimental and theoretical values of ot 
as a function of the K/tt ratio for the MINOS data. The er- 
rors in the experimental and theoretical values of ax are taken 
to be ± 0.012 and ± 0.013 respectively, obtained by com- 
bining the statistical errors in quadrature with the systematic 
errors in Tables and |III] but omitting the error in the K/tt 
ratio in each case. The error on the theoretical value of aT 
grows with increasing K/tt ratio because ex has a larger un- 
certainty than e^r, so a larger contribution from kaons intro- 
duces more uncertainty. The intersection of the two curves 
occurs at K/tt = 0.12tQ Qg. The uncertainty is estimated by 
assuming Gaussian errors for the the theoretical and exper- 
imental values of ar and performing a minimization to 
determine the Ax^ — 1 contour that encompasses the best fit 
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FIG. 7: The theoretical prediction for ar as a function of 
detector depth. The dashed (top) curve is the prediction using 
the pion-only model (of MACRO) and the dotted (bottom) 
curve is the prediction using a kaon-only model. The solid 
(middle) curve is the new prediction including both K and tt. 
These curves are illustrative only as the definition of effective 
temperature used to calculate the experimental values also 
depends on the K/tt ratio. The data from other experiments are 
shown for comparison only, and are from Barrett 1, 2 |2|1, 
AMANDA 1 4], MACRO [HJ, Torino Li2J, Sherman (Tj, 
Hobart Uj] andBaksan flTfl . 



the ± 0.014 TeV uncertainty in kaon critical energy; and e) 
the ± 0.003 TeV uncertainty in pion critical energy. These 
uncertainties are summarized in Table Hill 

Fig. I2] shows effective temperature coefficients from MI- 
NOS and other underground experiments, including those of 
the MACRO survey i3ll, as a function of detector depth. The 
MINOS and Sherman II15I1 effective temperature coefficients 
shown in Fig.|2]were calculated using Eq.[TOl The other exper- 



ulai' solid angle region on the sky and then normaUzing to the All-world 
Crouch underground muon intensity curve |26]. This was done with both 
Soudan 2 data |27] and with MINOS data |5], and these calculations were 
shown to agree to within 10%. Average cosmic ray muon flux, hke those 
determined here and in 1 5] can be determined using this method, although 
in any particular direction the rock map can be much different from what 
was calculated (e.g., in the direction of iron veins). 
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FIG. 8: The MINOS experimental qt as a function of the 
K/tt ratio (dot-dash line), with its error given by the 
cross-hatched region, on the same axes as the theoretical ax 
as a function of the K/tt ratio (dashed line), with its error given 
by the hatched region. The error on the experimental qt (from 
Table UllE'th cos 6, and (Tcff )) plus statistical error is 

± 0.012, and the theoretical ar error (from e,r,x and the rock 
map, Tablellllb is ± 0.013 at the best fit point. The 
intersection is at K/7r=0.12lQ gg. The solid line denotes the 
la contour around the best fit. 



point. 

Previous measurements of the K/tt ratio do not directly 
compare to this indirect measurement. Nevertheless, the cen- 
tral value of MINOS 's measurement is consistent with the 
collider-based direct measurements, although the indirect er- 
ror bars are larger than those on the direct measurements. A 
comparison of this measurement to other measurements is 
shown in Fig. |9l Only the MINOS result is for a reaction 
where the interacting particles do not have equivalent energy 
in the laboratory frame. Nevertheless, they are all presented 
on the same axes for a broad overview. The central value of 
MINOS' indirect cosmic ray-based K/ tt measurement is con- 
sistent with the collider-based direct measurements, and the 
associated error bars span the dispersion in those direct mea- 
surements. 



IV. CONCLUSIONS 

A five year sample of 67.32 million cosmic ray in- 
duced muons has been collected by the MINOS far de- 
tector and daily rate fluctuations have been compared to 
daily fluctuations in atmospheric temperature. These dis- 
tributions were shown to be highly correlated, with a cor- 
relation coefficient of 0.90. The constant of proportion- 
ality relating the two distributions, ar, was found to 
be 0.873 ± 0.009 (stat.) ± 0.010 (syst). This value 
is in good agreement with the theoretical expectation of 
(ar) — 0.864 ± 0.024. A measurement of the temper- 
ature dependence of the rate of /i+ separate from fi^ was per- 



FIG. 9: A compilation of selected measurements of K/tt for 
various center of mass energies. The STAR value was from 
Au+Au collisions at RHIC |29], the NA49 measurement was 
from Pb+Pb coUisions at SPS |30, 31], and the E735 
measurement was from p+p collisions at the Tevatron (33ll . 



formed for the first time. There is no statistically significant 
difference between Q!t(m^) ™d q;t(/^^). 

The experimental value of ar for the combined muon 
sample has the lowest uncertainty of any such measurement. 
While other experiments have estimated the effect of atmo- 
spheric temperature on kaon induced muons |12.|3|], this is the 
first result to quantify the kaon-inclusive effective temperature 
coefficient. The new kaon-inclusive model fits the MINOS far 
detector data better than the pion only model |3] and suggests 
a measurement of the atmospheric K/tt ratio. Applying the 
differing temperature variations of kaon and pion decay to the 
seasonal variations analysis allowed the first measurement of 
the atmospheric K/vr ratio for Ep >7 TeV. It was found to be 



K/tt = 0.12: 



-0.07 
-0.05- 
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